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Abstract 

Background: To obtain reliable quantitative real-time PCR data, normalization relative to stable housekeeping 
genes (HKGs) is required. However, in practice, expression levels of 'typical' housekeeping genes have been found 
to vary between tissues and under different experimental conditions. To date, validation studies of reference genes 
in pigs are relatively rare and have never been performed in porcine alveolar macrophages (AMs). In this study, 
expression stability of putative housekeeping genes were identified in the porcine AMs in response to the 
stimulation with two pathogen-associated molecular patterns (PAMPs) lipopolysaccharide (LPS) and lipoteichoic 
acid (LTA). Three different algorithms (geNorm, Normfinder and BestKeeper) were applied to assess the stability of 
HKGs. 

Results: The mRNA expression stability of nine commonly used reference genes {B2M, BLM, GAPDH, HPRTl, PPIA, 
RPL4, SDHA, TBP and YWHAZ) was determined by qRT-PCR in AMs that were stimulated by LPS and LTA in vitro. 
mRNA expression levels of all genes were found to be affected by the type of stimulation and duration of the 
stimulation {P < 0.0001). geNorm software revealed that SDIHA, B2M and RPL4 showed a high expression stability in 
the irrespective to the stimulation group, while SDIHA, YWIHAZ and RPL4 showed high stability in non-stimulated 
control group. In all cases, GAPDIH showed the least stability in geNorm. NormFinder revealed that SDIHA was the 
most stable gene in all the groups. Moreover, geNorm software suggested that the geometric mean of the three 
most stable genes would be the suitable combination for accurate normalization of gene expression study. 

Conclusions: There was discrepancy in the ranking order of reference genes obtained by different analysing 
algorithms. In conclusion, the geometric mean of the SDIHA, YWHAZ and RPL4 seemed to be the most appropriate 
combination of HKGs for accurate normalization of gene expression data in porcine AMs without knowing the type 
of bacterial pathogenic status of the animals. 
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Background 

Alveolar macrophages (AMs) are thought to be critical 
in the pathogenesis of several lung diseases [1]. Swine 
respiratory diseases, which has been described world- 
wide, affects swine of all ages and has a serious impact 
on economy, ecology and animal welfare in the pig rear- 
ing industry [2]. Both Gram-positive and Gram-negative 
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bacteria are causing respiratory disease in pigs [3]. As 
an in vitro model for the development of lung inflam- 
mation, AMs stimulation with PAMPs in culture is 
being frequently used for immunogenetic research in 
pigs [4-7]. Lipopolysaccharide (LPS) and lipoteichoic 
acid (LTA) are the PAMPs of the Gram-negative and 
the Gram-positive bacterial cell wall that cause activa- 
tion of an acute inflammatory response in vitro as well 
as in vivo. Gene expression assay is a common way to 
investigate the defensive role of AMs in the bacterial 
infections as well as to dissect the pathogenesis of 
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bacterial lung diseases. With this purposes, several stu- 
dies focusing on gene expressions have been conducted 
in AMs in vitro [4-7]. The gene expression are required 
to normalize for housekeeping genes (HKGs) which 
have tremendous effect on the results of expression 
study [8]. Therefore, it is crucial to know whether the 
expression stability of HKGs in AMs is affected by var- 
ious PAMPs from infectious agents but these data are 
currently unavailable for pigs. 

Quantitative real-time PGR (qRT-PGR) is a powerful 
technique for gene expression studies, which have 
become increasingly important in a large number of 
clinical and scientific fields [8,9]. Besides being a power- 
ful technique, inappropriate data normalization is the 
most important problem in qRT-PGR [8]. For an exact 
comparison of mRNA transcription in different samples 
or tissues, it is crucial to choose the appropriate refer- 
ence gene [9]. The most accepted approach to mRNA 
quantification is normalization of the expression level of 
a gene of interest (target gene) to the expression level of 
a stably expressed internal reference gene [8,9]. Normal- 
izing to a reference gene is a widely used method 
because it is simple in theory. Normalizing to a single 
reference gene is often used but Vandesompele et al. [8] 
suggested that geometric mean of multiple carefully 
selected HKGs is recommendable and suitable for accu- 
rate normalization. The normalization adjusts for differ- 
ences in the quality or quantity of template RNA or 
starting material and differences in RNA preparation 
and cDNA synthesis, since the reference gene is exposed 
to the same preparation steps as the gene of interest. 
This allows the direct comparison of normalized tran- 
script expression levels between samples. Reference 
genes should ideally be constitutively expressed by all 
cell types and should not be affected by disease and 
experimental procedure. To date, a universal reference 
gene has not been identified. HKGs are most commonly 
used reference genes [8]. Although HKGs are expressed 
by any cell, their expression varies among different cell 
types/organs, age, sex and treatment or experimental 
conditions [10-17]. Use of HKGs as reference genes for 
a particular sample type should be, therefore, validated. 

Ideally, the conditions of the experiment should not 
influence the expression of the reference genes [18]. 
However, the mRNA expression of reference genes from 
different cells and tissues [18-21] such as from AMs 
[1,10] may fluctuate due to infectious agents in vitro. 
Alveolar macrophages are being used as an important 
model to dissect the pathogenesis and genetics behind 
the infection through gene expression studies 
[5,6,22,23]. To date, no reference genes have been vali- 
dated for expression studies of AMs in pigs. The aim of 
this study was therefore to identify a set of stably 
expressed reference genes in porcine AMs cells 



irrespective of stimulation as well as in the case of sti- 
mulation by bacterial LTA and LPS in vitro. 

Methods 

Animals and preparation of alveolar macrophage cells 

Fourty-day-old three German Landrace piglets were 
euthanized for sampling. All animals were healthy and 
exhibited no signs of hypoxia or asphyxia or infections. 
Animals were kept and euthanized in the research sta- 
tion of Frankenforst at University of Bonn, following 
German pig breeding guidelines [24]. AMs were 
obtained from bronchoalveolar lavage (BAL) of animals. 
In brief, lungs were lavaged with 200 ml ice-cold sterile 
calcium-magnesium free Dulbecco's phosphate-buffered 
saline (PBS) (pH 7.4) that was instilled gently in 25 ml 
aliquots into the each of two adjacent lung subsegments 
and withdrawn immediately. BAL fluid from each ani- 
mal was collected in separate tubes and filtered through 
sterile gauze. Cells were centrifuged at 4°C for 10 min at 
400 X g. Pellets of bronchoalveolar cells were washed 
twice with sterile D-PBS at 250 x g for 10 min and 
resuspended in 2 mM L-glutamine-containing complete 
RPMI-1640 media (Sigma) supplemented with 10% fetal 
calf serum (Invitrogen) and containing antibiotics and 
antimycotics (penicillin, streptomycin and amphotericin, 
Invitrogen). The average purity of AM cells was 91% 
and other cells were mostly polymorphonuclear cells 
(8%) and remaining was lymphocytes. The cell viability 
was determined by Trypan blue dye exclusion method 
(> 98% in all cases). 

Stimulation of alveolar macrophage cells with LPS and 
LTA 

The cells were counted using Haemocytometer 
(AbCam) and concentration was adjusted. The AMs 
were plated in ultra-low attachment polystyrene 24- 
wells plate (CellStar) at 2 x 10^ cells in 1 ml medium in 
each well. Plates were incubated at 37°C with 5% CO2 
(Heraeus Instrument) for 48 h. After 1 h incubation, 
cells were stimulated with LPS of Escherichia coli 055: 
B5 (Sigma) (10 (ig per ml per well), LTA of Staphylococ- 
cus aureus (Sigma) (10 (ig per ml per well) and with 
both of LPS and LTA (10 (ig per ml per well). Cells 
were then collected at 1, 4, 8, 12, 24 and 48 h after sti- 
mulation for RNA extraction and stored at -80°C. For 
every time point non- stimulated control group was also 
included. 

RNA extraction and cDNA synthesis 

Harvested AM cells were washed in RPMI-1640 med- 
ium and the total RNA was extracted using Pico-Pure 
RNA isolation kit following the manufacturer's protocol 
(Arc turns. Applied Biosystems). In order to remove pos- 
sible contaminating genomic DNA, the extracted RNA 
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was treated with 5 [i\ RQl DNase buffer, 5 units DNase 
and 40 units of RNase inhibitor in a 40 [il reaction 
volume. The mixture was incubated at 37°C for 1 h fol- 
lowed by purification with the RNeasy Mini Kit (Qia- 
gen). Concentration of clean-up RNA was determined 
spectrophotometrically by using the NanoDrop ND- 
8000 (Thermo Scientific) instrument; the purity of RNA 
was estimated by the ratio A260/A280 with respect to 
contaminants that absorb in the UV. Additional exami- 
nation of integrity was done by denaturing agarose gel 
electrophoresis and ethidium bromide staining. Finally, 
the purified RNA was stored at -80°C for further analy- 
sis. Approximately 1.5 [ig of total RNA for each sample 
was transcribed into cDNA. cDNA was synthesised with 
SuperScript-II RT kit (Invitrogen). All samples were 
reverse transcribed under the same conditions. The 
synthesized cDNA was stored at -20°C and used in 
qRT-PCR reactions as a template. 

Selection of reference genes and primer design 

There are few previous studies for validation of selected 
HKGs across various tissues in pigs [11,12,16,25,26] 
with specific purpose and no study was devoted to vali- 
date reference genes in the AMs in case of inflammatory 
disease condition or in response to the bacterial product 
LPS and/or LTA. However, 'traditional' reference genes 
like GAPDH and TBP have been most often used in 
pigs [12,16,27-32]. Regarding porcine organs, ACTB, 
B2M, GAPDH, HMBS, HPRTl, RPL4, SDH A, TBP and 
YWHAZ have been previously compared [16,21]. More 



specifically in recent days, GAPDH, ACTB, RPL27, 
RPS29, RPS13 are compared in porcine stomach [31]; 
GAPDH, TBP, HPRT, RPS29, ACTB and RPL27 are vali- 
dated in porcine adipose tissues in different breeds of 
pigs [26] and B2M, SDHA, ACTB, GAPDH, HPRTl and 
TBP expression stability are compared in porcine mus- 
cle and liver tissues in pigs [25]. The genes used in our 
study were selected based on these previous studies. 
The following nine commonly used reference genes 
were selected: ACTB, GAPDH, HPRTl, B2M, SDHA, 
RPL4, YWHAZ, TBP and PPIA (Table 1). Primers were 
designed using the publicly available web-based Primer3 
program [33] and are listed in Table 1. They were tested 
using a BLAST analysis against the NCBI database 
http://www.ncbi.nlm.nih.gov/tools/primer-blast. 

Quantitative real-time PGR (qRT-PCR) 

Nine-fold serial dilution of plasmids DNA were pre- 
pared and used as template for the generation of the 
standard curve. In each run, the 96-well microtiter plate 
contained each cDNA sample, plasmid standards for the 
standard curves and no-template control. A no-template 
control (NTC) was included in each run for each gene 
to check for contamination. Quantitative real-time RT- 
PCR (qRT-PCR) was set up using 2 (il first-strand 
cDNA template, 7.4 \i\ deionized H2O, 0.3 (iM of 
upstream and downstream primers and 10 \i\ Ix Power 
SYBR Green I master mix with ROX as reference dye 
(Bio-Rad). The thermal cycling conditions were 3 min at 
95°C followed by 15 s at 95°C (40 cycles) and 1 min at 



Table 1 Selected candidate reference genes, primers, and PCR reactions efficiencies 



Gene 
name 


GenBank 
accession 
number 


Primer sequence 


Amplicon 
length (bp) 


Amplification 
efficiency (%) 


3^2 




Average Ct of cDNA 














Control LPS 


LTA 


Combined 


B2M 


NM_21 3978.1 


F:AC^CACACCGCTCCAGT 
R:CGGATGGAACCCAGATACAT 


180 


89.45 


0.992 


25.46 


24.30 


23.58 


23.34 


BLM 


NM_001 123084.1 


F:TCCTCACCTOTGCAmCC 
R:GTGGTGGCTGAGAATCCTGT 


152 


93.12 


0.993 


30.47 


28.58 


27.54 


28.06 


GAPDH 


AFOl 7079.1 


F:ACCCAGAAGACTGTGGATGG R: 
ACGCCTGCTOACCACCTO 


247 


89.45 


0.994 


36.96 


35.59 


33.92 


34.32 


HPRTl 


NM_001 032376.2 


F:AACCTOCmCCTOGTCA 
R:TCAAGGGCATAGCCTACCAC 


150 


91.88 


0.997 


29.21 


28.16 


28.50 


27.89 


PPIA 


NM_21 4353.1 


F:CACAAACGGI ICCCAGI 1 1 1 
R:TGTCCACAGTCAGCAATGGT 


171 


91.32 


0.997 


23.79 


22.88 


23.38 


23.27 


RPL4 


DQ845 176.1 


F:AGGAGGCTGTCTGCTOTG 
R:TCCAGGGATGmCTGAAGG 


185 


90.21 


0.993 


25.47 


24.92 


23.57 


24.19 


SDHA 


DQl 781 28.1 


F:AGAGCCTCAAGTOGGGAAG 
RGAGGAGATCCAAGGCAAAAT 


149 


92.24 


0.996 


30.35 


29.28 


28.21 


28.11 


TBP 


DQl 781 29.1 


F:ACGTOGGmAGGTOCAG 
R:GCAGCACAGTACGAGCAACT 


118 


99.43 


0.997 


31.81 


30.59 


31.11 


30.25 


YWHAZ 


DQl 781 30.1 


F:ATOGGTCTGGCCC™ACT 
R:GCGTGCTGTCmGTATGACTC 


146 


94.52 


0.994 


24.50 


23.74 


23.78 


23.23 



^ R^, correlation coefficient calculated from slope of the standard curve 
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60°C. Experiments were performed using the StepOne- 
Plus™ Real-Time PGR System (Applied Biosystems). 
Based on the Ct-values for all dilution points in a series, 
a standard curve was generated using linear regression 
and the PGR amplification efficiency of each primer pair 
is calculated from the slope of a standard curve [15]. 
Melting curve analysis was performed to verify the pre- 
sence of gene-specific peaks and the absence of primer 
dimmers (Figure lb). Agarose gel electrophoresis was 
performed to test for the specificity of the amplicons 
(Figure la). To ensure repeatability of the experiments, 
all the reactions were executed in duplicate and the 
mean was used for further analysis. 

Determination of reference gene expression stability 

The raw qRT-PCR amplification data was exported from 
the StepOne® software (Applied Biosystem) to Micro- 
soft® Excel. The averages of the Ct-values for each 
duplicate were used for stability comparison of candi- 
date reference genes in the NormFinder, geNorm and 
BestKeeper software. For easy understanding, the sam- 
ples were grouped into 5 different categories such as 
LPS stimulated, LTA stimulated, LPS + LTA (com- 
bined), control and irrespective to stimulation group 
(when all the stimulated and non-stimulated control 
were considered together). The effect of stimulation and 
time on the expression of housekeeping genes was 
tested using GLM procedure of the SAS software 
(ver.9.2; SAS, SAS Institute Inc., Gary, NG, USA). Differ- 
ences in gene expression levels between time and stimu- 
lation were determined using ^-test in SAS. P < 0.05 was 
considered statistically significant. 

Gt-values of all samples were exported to Excel, 
ordered for use in geNormPlus software (15 days free 
trial version qBasePlus; http://www.biogazelle.com) and 
geNorm transformed to relative quantities using the 
gene-specific PGR amplification efficiency [34]. These 
relative quantities were then exported to geNormPlus to 
analyze gene expression stability [8]. The approach of 
reference gene selection implemented in geNorm relies 
on the principle that the expression ratio of two ideal 
reference genes should be identical in all samples, inde- 
pendent of the treatment, condition, or tissue type. 
Increasing variations in the expression ratio between two 
genes correspond to lower expression stability across 
samples. geNorm calculates the stability using a pairwise 
comparison model [8]. geNorm determines the level of 
pairwise variation for each reference gene with all other 
reference genes as the standard deviation of the logarith- 
mically transformed expression ratios. In this way, the 
reference gene expression stability measure {M value) 
was calculated as the average pairwise variation of a par- 
ticular gene with all other control genes included in the 
analysis [8,15]. Lower M values represent higher 



expression stabilities. Sequential elimination of the least 
stable gene (highest M value) generates a ranking of 
genes according to their M values and results in the iden- 
tification of the genes with the most stable expression in 
the samples under analysis. geNorm was also used to 
estimate the normalization factor (NF„) using n multiple 
reference genes, by calculating the geometric mean of the 
expression levels of the n best reference genes [8]. The 
optimisation of the number of reference genes starts with 
the inclusion of the two genes with the lowest M value, 
and continues by sequentially adding genes with increas- 
ing values of M, Thus, geNorm calculates the pairwise 
variation VJ^n+i between two sequential normalization 
factors NF^ and NF„+i containing an increasing number 
of reference genes [8]. A large variation means that the 
added gene has a significant effect on the normalization 
and should preferably be included for calculation of a 
reliable normalization factor. Ideally, extra reference 
genes are included until the variation V„/V^+i drops 
below a given threshold. According to geNorm, if Vn/n+i 
< 0.15 the inclusion of an additional reference gene is not 
required and the recommended number of reference 
genes is given by n [8]. 

NormFinder uses an ANOVA-based model [35]. The 
software calculates a stability value for all candidate 
reference genes tested. The stability value is based on 
the combined estimate of intra- and inter-group expres- 
sion variations of the genes studied [35]. For each gene, 
the average Gt value of each duplicate reaction was con- 
verted to relative quantity data as described for geNorm, 
to calculate the stability value with NormFinder pro- 
gram [35]. The NormFinder reference tool was appUed 
to rank the candidate reference gene expression stability 
for all samples with no subgroup determination (irre- 
spective to stimulation) as well as with stimulation (LPS, 
LTA, and both LPS and LTA) as subgroup. A low stabi- 
lity value, indicating a low combined intra- and inter- 
group variation, indicates high expression stability [35]. 

The average Gt- value of each duplicate reaction was used 
(without conversion to relative quantity) in BestKeeper to 
analyze the stability value of studied genes [36]. BestKeeper 
creates a pairwise correlation coefficient between each gene 
and the BestKeeper index (BI). This index is the geometric 
mean of the Gt-values of all candidate reference genes 
grouped together. BestKeeper also calculates standard 
deviation (SD) of the Gt-values between the whole data set. 
The gene with the highest coefficient of correlation with 
the BI indicates the highest stability [36]. 

Results 

Purity, quantity of extracted RNA and verification of 
amplicons 

The optical density (OD) ratio A260/A280 nm measured 
with a Nanodrop spectrophotometer was 1.94 ± 0.17 
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M B21I BLM GAPDH HPRT1 PPIA RPU SDHA TBP YWHAZ 






Figure 1 Confirmation of amplicon size and primer specificity of studied genes, a) Agarose gel electrophoresis showing specific reverse 
transcription PCR products of the expected size for each gene, M represents DNA size marker, b) Melting curve analysis for all amplicons. 

V J 



(OD A260/A280 ratio ± SD). The average RNA concen- 
tration after extraction using Pico Pure was 10.33 (ig/(il 
±1.1 ((-ig/(il ± SD). The results of the averaged ampUfi- 
cation efficiencies are shown in Table 1. The amplifica- 
tion efficiencies for the nine candidate reference genes 
ranged between 89.45% and 99.43%. The agarose gel 
electrophoresis (Figure la) and melting curve analysis 
(Figure Ib-j and Table 1) revealed that all primer pairs 
amplified a single PCR product with expected size. 
Furthermore, sequence analysis of cloned amplicons 



revealed that all sequenced amplified fragments were 
identical to sequences used for primer design from Gen- 
Bank (data not shown). 

Expression levels of candidate reference genes 

Transcript abundance of commonly used HKGs were 
analysed in the different samples by direct comparison 
of their cycle threshold (Ct), assuming equal Ct for 
equal transcript number since all qRT-PCR reactions 
were performed with an equal quantity of total RNA. 
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Figure 2 showed that seven out of the nine selected 
genes presented Ct-values that ranged from 20 to 30 
cycles, while Ct-values from GAPDH (mean Ct 35.20) 
and TBP (mean Ct 30.94) were lower. The Ct of the 
remaining selected genes showed a reasonable disper- 
sion to moderately high expression levels. The expres- 
sion PPIA and YWHAZ were followed by B2M (mean 
Ct 24.17), RPL4 (mean Ct 24.54), HPRTl (mean Ct 
28.44) and BLM (mean Ct 28.67). GAPDH expression 
was lowest as indicated by Ct-values around 35 cycles, 
but it exhibited rather high dispersion over the stimula- 
tions and culture conditions indicated by large whiskers 
of the box (Figure 2a). According to variance analysis, 
the expression of eight genes was different from each 
other {P < 0.0001) (Figure 2). 

LPS and LTA affect expression level of reference genes 

The current study investigated fluctuations in expression 
of nine HKGs in AMs cultured with no stimulation, or 
stimulated with LPS, LTA or both. There were some 
fluctuations in the expression level of these genes in cer- 
tain conditions. The expression differences of these 
genes are shown in Figure 3. The variance analysis 
results between treatment groups and time of stimuli to 
the AMs are shown in Additional file 1: Table SI. Cell 
harvest time significantly affected the expression level of 
HKGs {P < 0.0001) (Additional file 1: Table SI). When 
the no stimulation control group was compared with 
the stimulated groups, the expression levels of all genes 
were lower in non-stimulated control group (Figure 3, 
Table 1). With LPS stimulation, mRNA expression levels 
of nine genes were increased. Beside in the case of LTA 
stimulation, expression levels of nine genes were 
increased compared to control group. Only expression 
of PPIA was decreased when cells were stimulated with 




Figure 2 Average cycle threshold (Ct) values of candidate 
reference genes tested in AMs under different conditions, a) 

The values are the average qRT-PCR cycle threshold numbers (Ct 
values). The bars indicate standard deviation. Letters indicate a 
significant difference in average Ct value. Average Ct values that 
have the same letter are not significantly different {P > 0.05) 

v J 



both LPS and LTA compared to LPS stimulation only 
(Figure 3e). 

identification of optimal reference genes 

Transcription profiling using qRT-PCR assays was then 
performed with these nine candidate genes, in samples 
from the four different conditions of AM cultures (LPS, 
LTA, combined LPS and LTA, and control). These raw 
Ct data were then analysed using different algorithms to 
identify the most suitable candidate genes. In each inde- 
pendent culture, the 9 genes were ranked according to 
their gene expression stability measure "M" (Figure 4a-e, 
left panel) with using the geNorm algorithm. Stepwise 
exclusion of the least stable gene allowed the genes to 
be ranked according to their M value (the lower the M 
value, the higher the gene's expression stability) [8]. All 
genes presented an M value below 1.5, which is the 
default limit for acceptable expression stability as 
defined by Vandesompele et al. [8]. Figure 4a shows the 
ranking of the nine candidate reference genes across the 
AMs based on their stability values without considering 
the type of stimulation of cells i.e. irrespective of stimu- 
lation group. SDHA, B2M and RPL4 were identified as 
the most stable HKGs (Figure 4a) in the irrespective of 
stimulation group. In case of the control group, geNorm 
showed that SDHA, B2M and RPL4 were the most 
stable HKGs (Figure 4b). When AMs were stimulated 
with Gram negative bacterial product LPS, geNorm 
identified B2M, SDHA and YWHAZ as the most stable 
HKGs (Figure 4c). YWHAZ, PPIA and RPL4 were the 
most stably expressed HKGs in the case of Gram-posi- 
tive bacterial product (LTA) stimulation group (Figure 
4d). When LPS was used combined with LTA for the 
stimulation of AMs, HPRTl, YWHAZ and SDHA 
remained the most stable genes (Figure 4e). All investi- 
gated groups identified GAPDH as the least stable refer- 
ence gene by geNorm (Figure 4a, c, d and 4e) except in 
control group where BLM was the least stable HKG 
(Figure 4b). 

NormFinder software ranked all HKGs according to 
their stability value (Figure 4f-j) [35]. The expression 
stability was not always consistent between the used 
softwares. By using NormFinder, genes SDHA, YWHAZ 
and HPRTl were ranked as the most stable HKGs in 
irrespective to stimulation group (Figure 4f). In the non- 
stimulated control group and LPS stimulated group, 
SDHA, YWHAZ and RPL4 remained the most stable 
genes (Figure 4g-h). In the LTA stimulated group, 
SDHA, YWHAZ and PPIA were ranked as the most 
stable HKGs (Figure 4i). In the combined LPS and LTA 
stimulation group, SDHA, HPRTl and TBP were found 
to be most stable HKGs (Figure 4j). PPIA remained the 
least stable HKGs followed by GAPDH and BLM 
according to the NormFinder algorithm. 
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NS LPS LTA LPS+LTA NS LPS LTA LPS+LTA NS LPS LTA LPS+LTA 



Figure 3 Expression levels of a representative subset of nine HKGs. a) B2M, b) BLM, c) GAPDH, d) HPRTl, e) PPIA, f) RPL4, g) SDHA, h) TBP 
and i) YWHAZ. Irrespective to stimulation: when all the stimulated and non-stimulated control were considered together; NS: no stimulation; LPS: 
lipopolysaccharide; LTA: lipoteichoic acid; LPS + LTA (combined): lipopolysaccharide used together with lipoteichoic acid. Differences among 
groups indicated with stars. *: P < 0.05; **: P < 0.01; ***: P < 0.001 



The results of reference gene evaluation by the Best- 
Keeper tool are shown in Table 2. According to the 
variability observed, candidate reference genes can be 
identified as the most stable genes, as they exhibited the 
lowest coefficient of variance (CV ± SD). A low SD of 
the cycle threshold (Ct) values should be expected for a 
useful reference gene. It is important to note that, genes 
that show a SD higher than 1 should be considered as 
unacceptable [36,37]. In the irrespective to stimulation 
group, YWHAZ was identified as the most stable HKG 
whereas GAPDH, BLM and B2M were removed from 
the initial statistics (Table 2). In the control group, 
PPIA was shown to be the most stable HKG by Best- 
Keeper. SDHA, BLM and RPL4 was identified as the 
most stable HKG by NormFinder in the case of LPS, 
LTA stimulated and combined LPS and LTA stimulated 
group, respectively (Table 2); whereas in all these cases, 
only GAPDH was eliminated from the initial statistics. 

Determination of the optimal number of reference genes 
for normalization 

The geNorm program calculates the normalization fac- 
tor assessing the optimal number of reference genes for 
generating the M factor by calculating the pair-wise var- 
iation K The pair-wise variation between these genes 



defines the variable V. The lower the variable V is, the 
less variation. The overall results are shown in Figure 5. 
For the irrespective to stimulation and combined LPS 
and LTA groups as shown in Figure 5a and 5e, five 
endogenous control genes are necessary to obtain the 
lowest changing V values in the analysed samples. On 
the other hand, seven endogenous HKGs were required 
for both the LPS and LTA stimulated groups (Figure 5c 
and 5d). For the control group, six HKGs were required 
to obtain an accurate normaUzation factor (Figure 5b). 
However, it is impractical to use excessive numbers of 
endogenous control genes for normalization, particularly 
when only a small number of target genes need to be 
studied or for rare samples that are very difficult to 
acquire [8,12]. Therefore, the use of the three most 
stable HKGs for the calculation of the NF was consid- 
ered acceptable for the majority of experiments [8,12]. 
To verify that the use of three HKGs simultaneously is 
adequate for normalization of qRT-PCR data, the corre- 
lation of NF values between the geometric means of the 
three most stable genes and the optimal number of 
genes was calculated for all sample groups. As shown in 
Figure 6, there is a high correlation between the two NF 
measures (i.e., the theoretical optimal number and pro- 
posed number, three) for all groups including 
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Figure 4 Ranking of nine candidate reference genes using geNorm and NormFinder softwares, (a-e) geNorm ranks the candidate 
reference genes based on their stability parameter M. The lower the M value, the higher the expression stability, (f-j) NormFinder ranks the 
genes based on a calculated stability value. The lower the stability value, the higher the expression stability. Irrespective to stimulation: when all 
the stimulated and non-stimulated control were considered together; Control: no stimulation; LPS: lipopolysaccharide; LTA: lipoteichoic acid; LPS 
+ LTA (combined): lipopolysaccharide used together with lipoteichoic acid. 
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Table 2 Expression stability of nine candidate reference genes evaluated by BestKeeper software 





B2M 


BLM 


GAPDH 


HPRTl 


PPIA 


RPL4 


SDHA 


TBP 


YWHAZ 


Bl 


Irrespective to stimulotion 


n 48 48 48 48 48 48 48 48 48 48 


SD [± Ct] 


1.05 


1.10 


1.61 


0.78 


0.79 


0.89 


1.03 


0.81 


0.76 


0.83 


CV [% Ct] 


4.36 


3.84 


4.57 


2.73 


3.37 


3.64 


3.57 


2.63 


3.20 


3.04 


Control 


n 


12 


12 


12 


12 


12 


12 


12 


12 


12 


12 


SD [± Ct] 


0.99 


1.33 


1.07 


0.93 


0.63 


1.07 


1.00 


0.74 


0.82 


0.93 


CV [% Ct] 


3.88 


4.35 


2.91 


3.19 


2.65 


4.18 


3.30 


2.34 


3.34 


3.27 


LPS 


n 


12 


12 


12 


12 


12 


12 


12 


12 


12 


12 


SD [± Ct] 


0.46 


0.73 


1.25 


0.61 


0.60 


0.75 


0.44 


0.73 


0.56 


0.59 


CV [% Ct] 


1.89 


3.52 


2.54 


2.16 


2.60 


3.00 


1.52 


2.39 


2.35 


2.16 


LTA 


n 


12 


12 


12 


12 


12 


12 


12 


12 


12 


12 


SD [± Ct] 


0.95 


0.42 


1.53 


0.61 


0.75 


0.73 


0.91 


0.70 


0.72 


0.74 


CV [% Ct] 


4.04 


4.51 


1.53 


2.15 


3.19 


3.11 


3.23 


2.27 


3.04 


2.75 


LTA + LPS 


n 


12 


12 


12 


12 


12 


12 


12 


12 


12 


12 


SD [± Ct] 


0.59 


0.39 


1.33 


0.46 


0.98 


0.36 


0.39 


0.49 


0.53 


0.34 


CV [% Ct] 


2.52 


1.37 


3.89 


1.66 


4.21 


1.49 


1.38 


1.62 


2.30 


1.26 



Descriptive statistics of nine candidate reference genes based on their cycle threshold (Ct) values. In the last column the BestKeeper (Bl) index is computed 
together with the same descriptive parameters for nine genes. Abbreviations: n number of observations, CV [%Ct] the coefficient of variance expressed as a 
percentage on the Ct level, SD [± Ct] the standard deviation of the Ct. Results from overall tissues irrespective of stimulation, non-stimulated control and different 
stimulations (LPS, LTA and LPS + LTA) are shown 



irrespective to stimulation group (r = 0.93-0.98, Pearson) 
(Figure 6a-e). This result demonstrates that the three 
most stable HKGs are sufficient for an accurate normali- 
zation of qRT-PCR data [8,12]. 

Discussion 

Using reference genes that have a stable expression 
between the compared groups is crucial in gene expres- 
sion studies. Several studies have shown that the use of 
different reference genes can change the outcome and 
conclusions of a study [13,19,38]. Ideally, the internal 
control gene for quantitative gene expression studies 
should not be influenced by the conditions of the experi- 
ment. However, our study showed that expression of the 
HKGs was affected by stimulation type as well as stimula- 
tion duration (Additional file 1: Table SI). Therefore it is 
generally recommended that the stability of HKGs is 
being validated prior to expression studies. There are 
some reports of the expression levels of HKGs in various 
cells and tissues and also of the methods used to analyse 
the stability of these genes. Recent research has demon- 
strated that the expression of HKGs may be altered due 
to state of the organ [21,39], age [17,21,26] and experi- 
mental conditions [18,20,40]. In the characterization of 
the course of an inflammatory reaction, quantitative real- 



time PGR has become a powerful tool for detection of 
inflammatory parameters, including cytokines and Toll- 
like receptors (TLRs). This tool is particularly useful in 
pigs since commercial species-specific antibodies directed 
against pig cytokines and TLRs are not commonly avail- 
able. To best of our knowledge, there has not yet been a 
detailed evaluation of HKGs in swine AMs. Moreover, 
there has not been a detailed study under different types 
of stimulation such as LPS, LTA and combined LPS and 
LTA that might be indicated Gram-negative, Gram-posi- 
tive bacterial infection or co-infection of both types of 
bacteria in vivo. Although no in depth studies are appar- 
ent in the AMs cells, there have been numerous research 
papers which have used single HKGs for normalisation of 
gene expression in AMs. These have included the use of 
HPRTl [5], GAPDH [41] and 18S rRNA [7] for normaU- 
sation of gene expression. As a consequence, in this 
study, we evaluated the gene expression stability of nine 
commonly used HKGs in porcine AMs, and furthermore, 
assessed their stability in states of different inflammatory 
models such as in response to LPS and LTA. 

In recent years, there have been a number of research 
papers and reviews evaluating the selection and effect of 
controls on normalised gene expression data in various 
pig tissues. Gu et al. [12] involved in the validation of 
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a) Irrespective of stimulation 
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^ 0.05 
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Figure 5 Determination of the optimal number of reference genes for normalization. The geNorm software calculates the normalization 
factor from an increasing number of genes (starting with at least two) for which the variable V defines the pairwise variation between two 
sequential normalization factors. The lower the pairwise variation, the better is the combination of genes for reference. V5/6 for example, shows 
the variation between the normalization factors of five genes in relation to six genes and shows that six genes is the combination providing the 
lowest pairwise variation. Irrespective to stimulation: when all the stimulated and non-stimulated control were considered together; Control: no 
stimulation; LPS: lipopolysaccharide; LTA: lipoteichoic acid; LPS + LTA (combined): lipopolysaccharide used together with lipoteichoic acid. 



20 common endogenous control genes in 56 fat- and 
muscle-type tissues. Nygard et al. [16] investigated a 
vast number of tissues for 10 HKGs. Studies focusing on 
more specific tissues, including the backfat, longissimus 
dorsi muscle [11], liver [25], adipose [26], stomach [31] 
and mesenchymal stem cells [42] are being reported in 
pigs. Taken together, it is very difficult to find a univer- 
sal' reference gene having stable expression in all cell 
types and tissues, and in particular to find reference 
genes that remain stable under different experimental or 
infectious conditions. According to the NCBI-PubMed 
statistics [12], GAPDH and ACTB are the two mostly 
used porcine HKGs. But they have been shown to vary 
considerably and are consequently unsuitable as refer- 
ence genes for normalization of gene expression analysis 
in many cases [43-45]. We appUed three software pro- 
grams to our data as complementary analyses to obtain 
the most suitable genes for our experiments. Both algo- 
rithms resulted in an overall comparable order of genes. 
Two of the three best genes were always presented by 
geNorm and NormFinder. Although BestKeeper [36] is 
found on the same principle as geNorm, not in every 
case both algorithms displayed overlapping suitable 
HKGs. 



In the present study, geNorm and NormFinder showed 
that SDHA, YWHAZ and RPL4 are the most stable three 
HKGs in the control (without any stimulation) group as 
well as in stimulation groups (Figure 4). Our results are in 
good agreement with Piorkowska et al. [26] who identified 
GAPDH and TBP as the least stable HKGs for the porcine 
adipose tissue. Beside, TBP was always found to be as a 
moderately stably expressed gene in this study. Nygard et 
al. [16] reported that RPL4, TBP and YWHAZ have the 
highest stability across tissues collected from healthy pigs 
which are somewhat consistent with the present study. 
Pierzchala et al. [25] recently reported that HPRTl and 
TBP are the most stable HKGs in porcine liver and in dif- 
ferent skeletal muscle tissues but it could be found that 
HPRTl and TBP is moderately stable through different 
experiments conditions in this study (Figure 4). Moreover, 
Svobodova et al. [46] estimated HPRTl has the highest 
stability while GAPDH was unstable across different por- 
cine tissues which are in agreement with our result for the 
GAPDH but not for HPRTl, Because HPRTl was found to 
be moderately stable in our experiment, except in com- 
bined LPS and LTA group. 

To our knowledge, there are only two studies evaluat- 
ing the stability of reference genes in AMs. One being 
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Figure 6 Correlation between the NF of most three stable and optimal number endogenous control. Pearson's correlations between the 
NFs of three endogenous control genes (NFS) and optimal number (six) of endogenous control genes (NFopt) for a) all samples irrespective of 
stimulation, b) non-stimulated control, c) LTA stimulated AMs, d) LPS stimulated AMs and e) LPS and LTA together used for stimulation. 
Irrespective to stimulation: when all the stimulated and non-stimulated control were considered together; Control: no stimulation; LPS: 
lipopolysaccharide; LTA: lipoteichoic acid; LPS + LTA (combined): lipopolysaccharide used together with lipoteichoic acid. 



in human AMs [1] and the other being in the horse 
[10]; no data is available on the stability of reference 
genes in AMs of other mammaUan species. Ishii et al. 
[1] reported that HPRTl is the most stable HKG, 
whereas TBP is the least stable HKG in both the LPS 
stimulated and non-stimulated AMs in human which is 
in good agreement with our result using geNorm. (Fig- 
ure 4b-c). Beekman et al. [10] used geNorm to validate 
the candidate HKGs and found that GAPDH, SDHA, 
HPRT and RPL32 were the most stably expressed genes 
in bronchoalveolar lavage cells of horses with inflamma- 
tory airway disease with corticosteroids treatment. In 
this study, SDHA was identified as suitable reference 
gene by using NormFinder through the experiments 
which is agreement with the report in horse [10]. 

According to the BestKeeper analysis software, in the 
irrespective to stimulation group YWHAZ was detected 
in accordance with the NormFinder and partially with 
the geNorm results (Table 2; Figure 4). SDHA was iden- 
tified as the most stable gene in both geNorm and 
NormFinder (Figure 4); however, BestKeeper identified 
this gene as unsuitable according to its algorithm cri- 
teria. In the control group, although PPIA was identified 



as a stably expressed HKG by BestKeeper (Table 2), this 
gene was identified moderately stable in geNorm and 
NormFinder (Figure 4). By using the three software 
algorithms similar results were obtained in LPS stimu- 
lated group, where SDHA was identified as the most 
stable HKG. In the LTA stimulated group, although 
BLM was identified as the most stable HKG by Best- 
Keeper, but showed very low expression stability in geN- 
orm and NormFinder. In case of the combined LPS and 
LTA stimulated group, RPL4 was found to be the most 
stable gene by BestKeeper (Table 2); however, this gene 
ranked as the fourth most stable HKG by geNorm (Fig- 
ure 4). Several studies previously reported similar discre- 
pancies for the findings of BestKeeper [15,31,37] and 
importantly, few studies followed the BestKeeper analy- 
sis method compared to geNorm and NormFinder. It is 
important to note that very similar discrepancies 
between the different algorithms have been observed in 
previous studies comparing statistical analysis methods 
[10,15,31,37,47]. 

However, we found that the first three most stable 
reference genes in most cases were consistent between 
the software geNorm and NormFinder, even if they 
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were not in the exact same ranking order. Similar find- 
ings are reported by previous studies in horse, human 
and plants [10,13,15,47]. Such discrepancy could be 
explained by genes' co-regulation. Indeed, co-regulated 
genes may become highly ranked independently of their 
expression stabilities with geNorm software [35]. More- 
over, NormFinder takes into account variation across 
subgroups, thus avoiding artificial selection of co-regu- 
lated genes by analyzing the expression stability of can- 
didate genes independently from each other [8]. 
However, no studies dealing with porcine reference 
genes stability used other analysis methods except geN- 
orm [11,12,16,26,31,42]. 

As described above, geNorm also provides a measure 
for the best number of reference genes that should be 
used for optimal normalization. In agreement with sev- 
eral previous studies, we postulate that the use of more 
than one reference gene allows for a more accurate nor- 
malization than the use of only one reference gene 
[8,12,35]. Based on a cut-off point for the lvalue, as 
described by Vandesompele et al. [8], a combination of 
the several most stable reference genes was calculated as 
being optimal for gene expression studies in control and 
PAMPs stimulated porcine AMs (Figure 5). However, as 
we described above and other studies [8,12] recom- 
mended, the combination of the most three stable genes 
are appropriate for accurate normalization. 

Conclusions 

In conclusion, this investigation found evidence that 
there can be variation in the expression of commonly 
used HKGs due to different PAMPs. Due to the new 
influx of data suggesting alterations in mRNA expres- 
sion according to bacteria type, we feel that beside ther- 
apy uses or experimental condition, there needs to be 
special consideration given to the selection of HKGs 
based upon the bacterial pathogen identification. This 
indicates that the choice of reference genes cannot be 
transposed from on study to the other without valida- 
tion for the specifics of each experimental protocol. 
Since different bacterial pathogens are cooperating in 
the respiratory tract as co-infection, our results will 
shed light on pathogenic or disease status of experi- 
ments. In general, we recommend using the geometric 
mean of SDHA, B2M and RPL4 to guarantee suitable 
normalization in across the AMs with unknown respira- 
tory pathogenic condition in pigs. Since in the most 
cases. Gram-negative and Gram-positive bacteria are 
observed together in respiratory diseases, HPRTl, 
YWHAZ and SDHA might be an appropriate set of 
reference genes for the gene expression normalization in 
AM studies. SDHA, YWHAZ and RPL4 could be sug- 
gested in case of AMs without any stimulation. This 
study offers an appropriate set of HKGs that might be 



used in the normalization of gene expression data in 
vitro cultured porcine AMs. 
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